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We describe the fringe packet tracking system used to equalise the optical 
path lengths at the Infrared Optical Telescope Array (IOTA) interferometer. 
The measurement of closure phases requires obtaining fringes on three 
baselines simultaneously. This is accomplished using an algorithm based on 
double Fourier interferometry for obtaining the wavelength-dependent phase 
of the fringes and a group delay tracking algorithm for determining the 
position of the fringe packet. The comparison between data acquired with 
and without the fringe packet tracker shows about a factor 3 reduction of the 
error on the closure-phase measurement. The fringe packet tracker has been 
able so far to track fringes of signal-to-noise as low as 1.8 for stars as faint 
as niH = 7.0. @ 2008 Optical Society of America 

OCIS codes: 120.0120, 120.2650, 120.3180. 



1. Introduction 

Observations performed with long-baseline ground-based optical/infrared interferometers 
are strongly affected by the turbulent atmosphere. Turbulence can reduce the visibility of 



1 



fringes in many ways as described by Porro et al^ for pupil-plane (or coaxial) beam combina- 
tion and Thureau^ for image-plane beam combination. Turbulence randomly modulates the 
phases of the fringes which can then become unusable for image reconstruction. Using three 
or more telescopes removes the atmospheric phase contamination. This is done through the 
closure-phase technique pioneered in radio astronomy^ and only recently applied to long- 
baseline optical interferometry.^ The necessary condition for obtaining meaningful closure- 
phases is that the three fringe packets must be detected within the coherence time. This is 
achieved by keeping the optical path difference (OPD) to a minimum. 

Prior to 2002, IOTA had only two telescopes and a single baseline, and the fringes were 
usually kept inside the scan interval manually by the observers. The installation of the third 
telescope at IOTA required an increase in the level of automation in the instrument, because 
manual tracking is not practical with three baselines to adjust. In particular the requirement 
to measure closure-phases necessitated a system capable of keeping the fringe packets in the 
centre of the scan using the existing hardware dedicated to acquiring data. Fringes must be 
acquired in the same coherence time in order to measure a closure-phase. This is especially 
important at IOTA where the bandwidths are frequently relatively large (15-25%) and the 
fringe packets quite narrow (3-6 fringes). 

We present here a simple and fast algorithm, which uses double Fourier interferometry 
(DFI)^ to extract wavelength dependent information from the fringe packet and calculate 
its group delay. The same algorithm was independently discovered by Tubbs^ and written 
as a final-year undergraduate project but never published. In this report Tubbs tested the 
algorithm through simulated data and data obtained from the COAST interferometer, but 
did not implement a working fringe-packet tracker. We only became aware of this work after 
our algorithm was routinely used at the IOTA. 

The time-domain interferograms, recorded by an infrared detector for the purpose of 
measuring the physical parameters of stellar objects, are used in order to track the position 
of fringes. A fringe packet is generated by pupil-plane interference of starlight from two 
telescopes, where the optical path difference around the zero-path-difference point is changed 
linearly with time. Each pair of the three telescopes generates two such interferograms, 
with complementary intensities. The goal is to control the optical paths such that these 
interferograms occur nearly simultaneously, so that the closure phase can be measured from 
them. To do this, the algorithm in the present paper is applied to each fringe packet, and the 
output of the algorithm is used to control the optical paths before the next interferogram is 
produced. Up to about 10 interferograms are generated per second, using triangular sweeps 
of the optical paths. No additional hardware is required. 

The application of the algorithm is not restricted to stellar interferometry but can be 
applied to all cases where broad-band interference fringes must be tracked. Other algorithms 
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for fringe packet tracking which have been tested at IOTA are the subjects of separate 
pubhcations7"^^ 

2. The Instrument 

IOTA is a long-baseline optical interferometer located at the Smithsonian Institution's Whip- 
ple Observatory on Mount Hopkins, AZ, comprising three 45-cm diameter telescopes which 
can be positioned at 17 stations on an L-shaped track, where the arms are 15 m toward 
the south-east and 35 m toward the north-east. IOTA operated with 2 telescopes from 
1995-2003, and 3 telescopes since February 2002. The interferometer,^^' has been used as 
a testbed for new cutting-edge technologies, ^^"^^ and has produced astronomy results in the 
3-telescope configuration. ^^"^^ 

The three beams arriving from the vacuum delay-line tank hit three dichroic mirrors which 
separate the visible and infrared light. The visible beam continues toward the star tracker 
servo system. In the implementation discussed here, the infrared beam is reflected toward 
three flat mirrors and then three off-axis parabolas which focus the three beams on three 
single-mode (H-band) fibers feeding the I0NIC-3T integrated-optics beam-combiner.^^ 

Interference is achieved inside the integrated optics component, resulting in three output 
pairs TT radians out of phase in intensity. The interference fringes are recorded while two of 
the dichroics are piezo-driven to scan a path of about +50 /im and -50 /im, respectively, 
in order to scan through the fringe packet in the three beams. The six combined beams 
are then focused on six separate pixels of the PICNIC array^^ and recorded as time series 
for science measurement. The same time series is used by the fringe-packet tracker. The 
path difference calculated by the packet-tracker is fed back to the piezo-scanning dichroics 
for a fast tracking response. The piezo scanners are off-loaded of their additional offsets 
every second, when a fraction of the error signal is sent to the short delay lines which are 
responsible for tracking the geometric delay caused by the rotation of the Earth. 

3. Calculating the Fringe Position 

3. A. Tracking the Fringe Packet Using Double Fourier Interferometry 

Our method of fringe packet tracking at IOTA calculates the group delay of fringes dispersed 
with DFI, which is used to obtain the wavelength dependent phase from the fringe packet. 
This is done by scanning the fringe packet over an interval greater than the packet length, 
where the spectral resolution is proportional to the scan length. The group delay tracking 
(GDT) method has been applied to interferometry since the very beginning of the field, 
when Michelson^^ used a prism for dispersing and acquiring fringes visually at the 20-ft 
interferometer. Labeyrie^'' used the same system and demonstrated fringe acquisition on a 
two-telescope interferometer. 
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Several systems have been proposed since then, for correcting the optical path.^^'^^ GDT 
(also called dispersed fringe tracking when applied to image plane interferometry) has been 
routinely used at several interferometric facilities. ^''"^^ In fact, at IOTA, GDT was selected 
as the original method of path difference monitoring in the visible, but the system was 
set aside in favor of making infrared observations. 

We will perform now a derivation of the algorithm for a two-telescope interferometer 
which combines the light in the pupil plane through a beam-splitter. With this setup the 
light intensity /(^) from the two complementary outputs of the beam-splitter is measured 
by square-law detectors: 

where It is the mean intensity which becomes split by the beam combiner (± is for com- 
plementary outputs) and is modulated by the interferometric signal, function of the optical 
path ^. V is the contrast of the fringes. A/ (cm~^) the bandwidth of the ideal (rectangu- 
lar) spectral filter, /o (cm~^) the frequency of the fringe, (rad) a generic phase, and 
(cm) is the fringe-packet centre, which varies from scan to scan owing to atmospheric path 
fluctuation. This is the actual value we want to calculate in order to correct the optical path. 

In practice we are dealing with discrete intensities, expressed as a finite series of data, so 
we can rewrite (1) as the discrete function n{j), which is the detector count. 

"^^■^ = "° { ' ^ ^ '"Ja-Jotit""^ - + } + (2) 

Here uq is the mean number of electrons per channel, and V is the fringe visibility. The sample 
number j ranges from 1 to A^, where is the total number of samples; typically A^ = 256. 
The centre of the fringe packet is at sample number Jq. If the high- and low-wavenumber 
limits of the spectral filter are written as ruh and mi (waves per sample), then the filter 
full-width at half-maximum is Am = rrih — mi, and the filter centre is mo = (m/j + m;)/2, 
both in units of waves (or fringes) per sample. The phase (p (radians) is an offset between 
the sinusoid carrier wave and the centre of the sine envelope. The ej term represents additive 
noise, the sum of electron counting statistics, detector read noise, and random atmospheric 
scintillation. The goal is to extract an estimate of Jq from the ensemble of A^ data points. 

The first step is to calculate the fast Fourier transform of the data string n{j), and scaling 
the result so as to eliminate uninteresting factors. The result is a sequence of A^/2 complex 
numbers which can be written as 

h (mo) = exp (-i27rmo Jq) + rjmo (3) 

where 77^0 is the scaled transform of the noise term ej. 
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The second step is to select two of the n values, say n{mQ) and n(mo + Ami2), and form 
the cross-spectrum product of the first with the complex conjugate of the second, i.e., 

X (mo, mo + Ami2) = n(mo) ■ n (mo + Ami2)* . (4) 

As can be seen from (3), the cross-spectrum X will generate a complex number (plus noise) 
whose phase contains the unknown quantity Jq, multiplied by known terms. The additive 
noise would degrade our estimate of Jq from the cross spectrum term, however adding several 
such terms will improve the signal-to-noise ratio of the resulting complex term,^^ and likewise 
improve our estimate of Jq. 

The third step is to calculate the average cross-spectrum X, as for example in 

— ^ X(mo,mo + Amia) ^ . j . . _ 

X= —— __-_ = exp(-z27rA™i2Jo)+r7, (5) 

mo=mi y^h -mi + i) 

where rj is the new noise term, presumably smaller than the root-mean-square of r^^o by a 
factor on the order of ^ {rrih — mi^i). The limits on the sum are only suggested values, and 
could be changed so as to improve the signal-to-noise ratio if applicable. 

(We note that the cross-spectrum is also used in the Knox-Thompson algorithm^^ for 
recovering near diffraction-limited images of stellar objects from snapshots degraded by at- 
mospheric seeing.) 

The estimated fringe packet position < Jq > can now be recovered from the average 
cross-spectrum using 

arg (It. 



° 27rAmi2 

Here the estimated packet center < Jo > is measured in units of sample numbers from the 
first point in the interferogram. In practice we subtract N/2 from this value, and send the 
resulting value (suitably scaled) as an error signal to the servo system, such that the starting 
point of the next scan is adjusted accordingly. 

We note that any value of Ami2 could be used, but we suggest that Arriu = 1 is optimum 
in the sense that it is least likely to produce a complex S oi S that is biased by wrap-around 
effects, which will occur if the phase shift between adjacent values of n are separated by 
more than vr. Hereafter we drop the expectation value brackets, and write Jo for < Jq > for 
notational simplicity. 

3.B. Algorithm summary 

1. The interference fringes are recorded while two of the optical paths are piezo-driven to 
scan a path of about -|-50 fim and -50 fim, respectively. The time series nj is recorded 
by the infrared camera for the three pairs of beams. 
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2. The fast Fourier transform n{mQ) of the time series nj is computed separately for the 
three pairs. 

3. The cross-spectrum X is calculated for the 3 beams by multiplying ri(mo) by its 
complex conjugated n*{mo + 1) shifted by one sample. 

4. The average cross-spectrum X is computed in correspondence of the fringe peak in the 
cross-power-spectrum, in order to average complex vector with higher signal-to-noise. 

5. The position of the fringe packets is obtained from the phase of the average vector X. 
using Eqn. 6 

3.C. Baseline Bootstrapping 

With baseline bootstrapping^^' we are capable of blind-tracking fringes on a baseline when 
the signal-to-noise of the fringes is too low, provided that good signal-to-noise is available on 
the other two baselines. For this reason, we calculate the optical path Jq for three baselines 
even if we correct the path for two baselines. We can express one optical path as the weighted 
average of the other two optical paths, the weights being equal to the SNR for the fringes 
obtained on those baselines. The signal-to-noise is calculated from the cross spectrum \X\ , 
dividing the averaged power inside the fringe peak by the averaged power outside the fringe 
peak. We then observe that the optical path, in a closed loop must be equal to zero: 



Jol + Jo2 + Jo3 = (7) 

where JqI and Jo2 are the optical path where the servo loop is acting, while JqS is the 
reference optical path. To the path JqI, Jo2 and JqS are associated the weights wl, w2 and 
w3 respectively. We have two values for each optical path. One is the value obtained directly 
on that baseline (for example JqI with weight wl), the other is the value calculated from the 
linear combination of the other two baselines (for example JqI' = —Jo2 — Jq3 with weight 
wl' = {w2w3)/{w2 + w3)). The weighted average of Jq and Jq can then be written as:^^ 

— wlJol + wl'Jol' 
wl + wV 

where JqI is the weighted-averaged path difference. Similarly for Jq2: 

'Y^ — w2.]o2+w2' Jo2' 
"'O^ — w2+w2' 

Jo2' = -Jol-Jo3 (9) 



w2' 



wlw'i 



uil+u;3 

The advantage of using a weighting system is that we do not have to select the best 
baseline values a priori, but rather the weighting allows them to be selected automatically. 
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4. Simulations 



The linearity of the fringe packet tracker (FPT) algorithm was tested through simulations. 
This simulation was performed in presence of photon and detector noise, with an average 
number of 300 photoelectron/sample for the fringe intensity (Poisson distributed), a fringe 
visibility of one and an additional 12-electron mean of Gaussian noise for the detector. 

It was found that the algorithm is capable of detecting the sign of the correction, necessary 
to bring back the fringes to the centre of the scan even if the main part of the fringe packet 
is outside the current scan. The necessary condition is to have enough signal-to noise so that 
the side lobes are detectable. This behavior is shown in the graph of Fig. 1. 
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Fig. 1. The optical path returned by the algorithm with respect to the fringe 
position (continuous-line plot). The two vertical dotted-lines represent the in- 
terval inside of which the fringes are sampled, for a scan length of 50 /im. 
Three fringe packets (dashed-lines) mark three representative positions on the 
plot: outside the scan range (top), at the centre of the scan (centre) and half 
outside the scan range (bottom). Note that only the side lobes are visible in 
the scan range of the plot at the top. 



7 



The returned position depends linearly on the fringe positions when the packet is inside 
the scan but it is non-linear when the packet is outside the scan. Nevertheless the information 
of the side lobes can still be used to bring the fringes back to the centre of the scan. The 
algorithm is sensitive to the slope of the phase across the bandwidth and the sign of the 
correction is preserved. The algorithm starts to fail at the limit of the range where the signal 
is too low and the side lobes not detectable for the given integration time. The returned 
position depends linearly on the fringe positions when the packet is inside the scan and 
becomes strongly non-linear (but with the correct sign) when the packet is outside the scan. 

, O I I I 1 .0 1 — ' ' ' ' I 




Fig. 2. Closure phase measurement for the star WR140 {rriH = 5.3). The 
points represented on the plot are complex vectors normalised by the signal- 
to-noise. Each point is calculated using a single fringe measurement for each 
baseline (notably the N-S, N-W and W-S baselines, where the telescopes 
where positioned at N=35m, S=15m and W=10m). There are 200 points in 
the diagram. The vector represented in solid line is the average closure phase 
of the previous 200 vectors and the dashed line represents the error on the 
closure phase. Finally the left panel represents a closure phase measurement 
in open-loop mode (FPT not active, closure phase = 169.0 ± 4.3 deg) while 
the right panel is the closed-loop case (FPT active, closure phase = 168. 4± 1.5 
deg). 



5. Results 

The FPT has been routinely in use at the IOTA interferometer since early 2002. Fig. 3 shows 
the reduction of tracking residuals operated by the algorithm. A reduction of about a factor 
2 in the RMS change of the optical path, between data acquired with tracking switched off 
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compared to data with tracking on, can be shown from the recorded data. This hes mostly 
in the low frequency zone because of the limited bandwidth of the combined fringe- sensor, 
fringe-actuator. For the IOTA interferometer this is a significant improvement, since the high 
frequency phase noise is dealt within post-processing using the closure-phase information. 
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Fig. 3. The change of the OPD with the FPT turned off (top) and on (bottom). 
This experiment was performed using a single, 21 m long baseline at 1.65 /xm, 
during average seeing and no baseline bootstrapping. The data was recorded 
the 6th of March 2002. We measure a factor 2 reduction in the RMS change 
in optical path. 



It is in fact more important to be able to maintain the fringe packet superposed and being 
able to do so for faint sources rather than reduce the residuals of the OPD to a smaller value. 
Fig. 2 shows two measurements of closure phase. When the FPT is switched off the closure 
phase is 169.0 ± 4.3 deg but it is 168.4 ± 1.5 deg when the FPT is operating (a factor 3 error 
reduction for the closed-loop case). 

We also show that the FPT is capable to track fringes on a 7.0 mn star with a SNR as 
low as 1.8 as shown in Fig. 4. 
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Fig. 4. Faint fringe observed on a 7.0 mn star, yielding a SNR of 1.8. The 
dashed hne across the fringe shows the expected position calculated by the 
FTP. 

6. Conclusions 

We describe the fringe packet tracking system now routinely in use at the IOTA interferome- 
ter. The packet tracker uses existing hardware to perform its functions, notably the infrared 
camera used for data acquisition, as a fringe sensor and the fringe scanning platforms com- 
bined with the delay lines, for correcting the optical path. It is based on an algorithm which 
exploits the wavelength-dependent information contained in the fringe packet, which is ex- 
tracted using double Fourier interferometry. The phase of the fringes obtained for different 
frequencies is used for calculating the group delay of the fringe packet. Since we use baseline 
bootstrapping we are capable of blind-tracking fringes on a baseline when their signal-to- 
noise is very low, provided that the fringes have a good signal-to-noise on the other two 
baselines of the interferometer. 

We use numerical simulation to model the case of a fringe packet outside the scanning 
range, operated by the algorithm, when only the side lobes of the fringe packet are visible. 



10 



We also show that the fringe packet tracker algorithm effectively reduces the slowly varying 
atmospheric and instrumental-induced additional path, leaving the fast phase noise to be 
dealt with post-processing. Moreover the FPT delivers about a factor 3 reduction of the 
error on the closure-phase measurement and has been able so far to track fringes with SNR 
as low as 1.8. 
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